Superfluid transition in a rotating resonantly-interacting Fermi gas 
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We study a rotating atomic Fermi gas near a narrow s-wave Feshbach resonance in a uniaxial harmonic trap 
with frequencies Oj_, fl z . Our primary prediction is the upper-critical angular velocity, u C 2(S,T), as a func- 
tion of temperature T and resonance detuning 8, ranging across the BEC-BCS crossover. The rotation-driven 
suppression of superfluidity at ui C 2 is quite distinct in the BCS and BEC regimes, with the former controlled by 
Cooper-pair depairing and the latter by the dilution of bosonic molecules. At low T and Q z <C Q±, in the BCS 
and crossover regimes of < 8 < S c , w C 2 is implicitly given by h\J lj 2 2 + Q 2 ^ ~ 2A^/hftj_/ep, vanishing as 
lu C 2 ~ S1_l(1 — 8/Sc) 1 / 2 near 8 C « 2e F + ln(e F /Kl±) (with A the BCS gap and 7 the resonance width), 
and extending the bulk result hcj C 2 ~ 2A 2 /e F to a finite number of atoms in a trap. In the BEC regime of 
5 < we find uj C 2 —* fij, where molecular superfluidity can only be destroyed by large quantum fluctuations 
associated with comparable boson and vortex densities. 



PACS numbers: 03.75.Fi, 03.75.Lm. 05.30.Jp 

Recent advances in atomic gases near a Feshbach resonance 
(FR) have led to the realization of resonantly-paired atomic 
superfluids [1, 2]. The proximity to a FR allows a tunability of 
the pairing interaction, thereby permitting unprecedented ac- 
cess to fermionic superfluidity ranging from a weakly-paired 
BCS regime to a strongly-paired molecular BEC regime. 

A fundamental aspect of a superfluid is its non-classical 
response to an imposed rotation. Unable to exhibit rigid 
body rotation, a superfluid rotates by nucleating a vortex ar- 
ray with density n v = mw/Trh set by the rotation rate u>. 
Although considerable progress has been made in elucidat- 
ing the properties of rotating superfluids in the BEC regime 
(with bosonic atoms) [[3], 3, @], considerably less is under- 
stood for a resonantly-paired trapped superfluid in the BCS 
and crossover regimes. Recent spectacular experiments by 
Zwierlein, et al. [6] provide strong motivation for a study of 
these regimes. 

Additional motivation is provided by the relation of a rotat- 
ing superfluid to a type-II superconductor in a magnetic field 
iHM @], with the Coriolis force 2mr x u> in the former corre- 
sponding to the Lorentz force in the latter, with the identifica- 
tion of — eB I c with 2mu>. Although significant insight can be 
gained from this connection, it is limited to the BCS regime 
and does not include important ingredients that are unique to 
a rotating superfluid. These include: the absence of screen- 
ing, the centrifugal force, the trap potential, the fixed number 
of atoms, and the tunable resonant pairing interaction, all of 
which are absent in the analogous superconductor problem. 
In the latter a concomitant Zeeman field also appears. It can 
be effectively introduced into the atomic problem (but we will 
not do so here) by imposing a difference in the number of the 
two pairing atomic (hyperfine state) species [j iol] . 

In this Letter we study the effect of an imposed rotation lu 
on a trapped (with axially symmetric trap frequencies f2j_ and 
fi z and n = (f^Qj.) 1 / 3 ) resonantly-paired superfluid near 
a narrow FR (i.e., the width of the resonance, V, is smaller 
than the Fermi energy ep = (3N) l / 3 ffll = h 2 k F /m, 7 = 
^/T/eF < 1), tuned through the BEC-BCS crossover|fl5| . 



Our primary prediction is the upper-critical angular velocity 
Wc2(8, T) as a function of FR detuning S and temperature T, 
illustrated for T — > in Fig. Q] and for a range of (5 as a func- 
tion of Tin Fig. 
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FIG. 1 : (Color online) The upper critical rotation rate uj C 2 (8,T — > 0) 
as a function of the detuning 8 in an isotropic trap with 10 7 atoms, 
and for two FR widths 7 = 0, 1/3. uj C 2 vanishes for 8 greater than 
5 C w 2e F + 2 e F In(e F /ftfi±). 



For low T and Sl z <C in the BCS and crossover 

regimes of < 8 < 8 C we find u> C 2 to be implicitly given by 
hyj oj1 2 + Q,^_ fts 2A UJc2 ^hQ±/e F c2 , controlled by the dis- 
creteness of the rotated trap spectrum cutting off the weak 
Fermi-surface Cooper pairing characterized by the BCS gap 
A^o = A(c_f) = 8e~ 2 e F exp[- ^j== }■ In this regime, for 
I ln(e F /M]jJ > 1, this leads to 



,T=o 

J c2 




(B 6 



for 5 — ► S c , 
for < 6 < S c , 



(1) 



f i- 5re_F ln(eF/fi^±) cor- 




where the behavior near S c 2e 

responds to a vanishing of T" =u when the BCS conden 
sation energy A 2 /ef drops down to the trap level spacing 



2 



Ml±. At low detuning S <C S c , oj C 2 rises up to but is lim- 
ited below Qj_ by the implicit dependence entering through 
€ f = e F (l - uj 2 /fi^) 1 / 3 and A w = A(e F ) due to the cen- 
trifugal and Coriolis forces reducing the effective trap poten- 
tial and atom density. In this limit, we recover the bulk result 
Tuu C 2 rs 2A^ 2 /e^ c2 , corresponding to Gorkov's (fixed chem- 
ical potential) prediction for type-II superconductors [81]. 

In the opposite limit of 3 hi{e F "C I, S c ~ 2e F and 
w c2 =0 reduces to 



w<r 2 =0 (<5) ~ w *0*) = n ±y/l - (<V 2 <^) 3 , for < <S < 2e F , 

(2) 

This also smoothly matches onto the ui C 2 — > f2j_ behavior 
in the BEC regime of 5 < 0, where molecular superfluid- 
ity can only be destroyed by large quantum fluctuations as- 
sociated with comparable boson and vortex densities, where 
it undergoes transitions to a variety of bosonic quantum Hall 
states | llll . 

As can be seen in Fig. [J and Eq. (Q, a distinction between 
BCS (S ^> 2e F ) and crossover (0 < 5 < 2e F ) regimes is 
not reflected in lu C 2 (S), controlled throughout by weak Cooper 
pairing. However, T c (uS) (Fig. |2) does distinguish between 
these regimes, exhibiting a point of inflection at a scale ($) 
in the crossover (but not in the BCS) regime. At a rotation 
rate ui < uj*(S) (nonzero only in the crossover regime) a fi- 
nite number of molecular bosons N B (w) is present and T c (oj) 
is set by T® EC [N B (u), S)]. At a higher rate, w*(S) < uj < 
uj C 2(S) (as a consequence of an increased atomic density of 
states) e F drops below 5/2 and T c (ui) is determined by an ex- 
ponentially reduced bulk BCS value T c BCS (e F ) = ^A(e F ) 
(c « 0.577 is the Euler-Mascheroni constant). 
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FIG. 2: Main: T c (5, u>) in a spherical trap for a range of rotation 
rates u as a function of S. Inset: T c (5,u>) for a range of detuning 
S = {0, 0.9, 1.8, 2.7}ef as a function of u> with the top line being 
for S = 0. In both cases, JV = 10 7 and 7 = 1/3. 



In the BCS regime (S > 2e F ), at slow rotation u> <C uj^ , 
we find 



T c {u) _ 7C(3) ej-fil 



T^BCS 
^c.bulk 



127T 2 



(3) 



with the trap responsible for the quadratic in ui suppression of 
T c in the BCS regime, contrasting to the linear \ui\ suppres- 
sion in a superconductor [7]. In the BEC regime, the rotation 
suppresses T c according to 



T c (w) 



T BEC / 1 



2 /nl) lf \ for 6 <0. (4) 



We now sketch the derivation of these results, delaying details 
to a future publication jl2\\ . 

A gas composed of two hyperfine species (labeled by a =| 
, |) of fermionic atoms interacting through a tunable (via a 
"bare" detuning So) s-wave FR, corresponding to a closed- 
channel diatomic molecular state, is welFcharacterized by a 
Hamiltonian H = U F + U B + U g , [0, [H, 0] where 



h 2 V 2 
2m 



V(r)-fi-ujL z 



H B = fd 3 v^(v) 

H g =g f d 3 r(0t(r)^ T (r)^(r)+^](r)^}(r)^(r)) , (5) 



ft 2 V 2 
Am 



+ 2V(r) ~2fi + 5 -ujL 



Vv(r), 



written in the frame rotating with uiz. Here V'J( r ) is an atomic 
creation operator in the hyperfine state \o) whereas <fi'(r) is 
the molecular creation operator. Above, m is the fermion 
mass, fi the chemical potential, g the atom-molecule intercon- 
version amplitude, L z = —ihz-(v x V r ) the angular momen- 
tum operator, and V(r) — \m (fi^x 2 + y 2 ) + Q 2 z 2 } is the 
atomic trap potential. 

To determine uj C 2 (S, T) we look for an instability of the 
normal state to a superfluid state. To this end we consider 
system's partition function Z = Tr[e~" n ], (j3 = 1/fcsT) 
and integrate out Fermi fields perturbatively in the FR cou- 
pling g, valid when the corresponding dimensionless cou- 
pling 7 = ^^(^f 2 « 1. We obtain Z = 

J D^D0exp(-S , 0/S) with (defining h B = -h 2 V 2 /Am + 
2V{r) -2^ + S ~ uL z ), 



S,h 



•pn 



AO) 



drj d 3 r^(r,T) (hd T + h B ) 4>{t,- 
1 f pn f - 

-J drdr' J d 3 rd 3 r' <j>(r, r)E(r, r', r - r')cj){r' , r'), 



(6) 



- 5 2 ^ UJ (r,r',r)^ 0) (r,r',r) is the 



where E(r,r',r) = ~g 2 g ( F ] ' 
molecular (Cooper pair) self-energy arising from molecule 
fluctuations into a pair of atoms governed by the free 
fermion propagator <^ 0) (r, r', t) = -(Vv(r, r)Vv(r', 0)) 
and zf) — Tr [e~' 37iF ] is the partition function of free 
fermions. Higher-order corrections in g lead to the inclusion 
of molecular interactions induced by the atomic FR, small for 
7 < 1. 

The instability to the superfluid state at oj C 2 (S, T) is deter- 
mined by the vanishing of the lowest eigenvalue of the inverse 



3 



of the effective molecular propagator Q B 1 {r, r'; iu>i)\i ue =o = 
6^ 0)_1 (r, r'; i^)-E(r, r', i^)| ia)<= o, where (r, r', r) = 
— (0(r, r)0(r', 0))o is the free molecular propagator and oj£ = 
2£tt/(/3H), (I e Z) is the bosonic Matsubara frequency [16]. 
To lowest order in g 2 , ^ i3 1 (r, r'; 0) is diagonalized by eigen- 
states of 1 (r, r';0) (valid for a narrow resonance), 

that are simply rotated trap eigenstates n (r) defined by 
(-fr 2 V 2 /4m + 2V(r) - wL z )^„(r) = e°</>„(r), with eigen- 
values (n = {n , n, n z }) 

= htt ± (2n + 1) - 2ftcjn z + M7 Z (n + 1/2) , (7) 

where no £ N,n £ and n z = {— n, — n + 1, . . . , n} 



are the axial, radial and angular momentum quantum num- 
bers. The single-particle fermionic and bosonic spectra are 

t , • , 0{F) o j 0(B) i x o 

then given by e„ = e„ — /i and e n = e„ + oo — 2^i, 
respectively. In the weakly-interacting limit the condensation 
is therefore into n = = {0,0,0}, a nonrotating bosonic 
eigenmode localized at the center of the trap [HO. 



Projecting the self-energy operator S(r, r'; 0) onto the con- 
densate 0o (r) gives oj C 2 as the solution of the Thouless crite- 
rion e { B} = f° (B) - J d 3 rd 3 r'0 o (r)£(r, r'; O)^o(r'), explic- 
itly given by 



g , j tanh( /3e r ) /2) + tanh(/ 3e ^)/2) , , 3fc 

-Kn.n' ;r— — nTTvTx / TTT^TTT^ } > (°) 







where we have expressed the "bare" detuning, So, appearing in the Hamiltonian in terms of the physical detuning d = 5i 
g 2 J prp- arising in the two-body s-wave scattering measurement through the scattering length a s — —j^§rg- The matrix 
element R n . n > = 2~ 3 / 4 / d 3 r^ o (r)0* (r/\/2)0*/(r/\/2) is given by 

\ R | _ 1 (mn\ m [l + (-l)"°+"o](n + n )! (n + nQ! <? n , +n ,, 

n ' n '' 2 2 ("^)^)!y^! 2"+"' _ n z )l(n + n z )l(n> + n> z )l(n> ~ <)\ ' 



r 



Equation (8) needs to be solved together with the 
total atom number equation N = (3~ 1 d\nZ/dfi — 
J d 3 r (2{^(p) + (t/iJt/v)), fixing /i in terms of N, which to 
this order in g is given by TV = N F + 2N B + 0(g 2 ) with 



4 = 2E n #« 



+ l)and7V B = E n l/(f 



Analytic analysis of Eq. ® is possible due to considerable 
simplifications in the limits of a narrow resonance (7 <C 1) 
and a Fermi energy that is large compared to the rotating 
trap level spacing, tp 3> + w). The latter ensures 

that the fermionic trap level sums are dominated by large 
quantum numbers n, n', allowing the Gaussian approximation 
(§)«+»' (n+n')!/(n!n'l) « (TrnJ-^exp [-(n - n') 2 /4n] . 
Physically this approximation corresponds to a weak Coriolis 
force, with the atomic trajectories nearly straight lines (locally 
well-characterized by plane waves) turning by a small fraction 
of the particle spacing, n -1 / 3 . This reduces the Thouless cri- 
terion, Eq. dD, to 



F{fln) + G 



(10) 



where 

G(x, a, b) 



oc />oo />1 

/ dydp / dz 

!> JO JO 1 



p 2 

V 7TX 2 ' 



y-x 



cosh 2 (^) 
sinh 2 (f ^[a+(6-a)2 2 ]J 




/■°° -e~ p2 f 1 




/ dp — -=— I dz 


Cl + "00 Q + «<2^ 


/-oo V 71 " Jo 





tanh((y - x)/2) 1 



(11) 
(12) 
(13) 



y-x 

F(x > 1) = ln(8a;/7r) - 2 + c, 

with Q = py/x(a + (b — a)z 2 )/(27r), ci = In 4 + c, and 
ipo(y) the digamma function. Within the same approximation, 
the number equation reduces to 

3 



N 



-iis (- 



Li 3 



(14) 



where Lis(x) — YlkLi fr is the trilogarithm function, 
with asymptotic forms Li^{—e z ) w — z 3 /6 for z 3> 1, 
Li 3 (-e z ) « -e 2 for z < -1, and £i 3 (l) = C(3) « 1.202 . 
The first and second terms are the number of atoms and ther- 
mally excited molecules (with the condensate vanishing at 



4 



oj C 2), respectively, and = Q(l — o^/fijj 1 ' 3 is the ef- 

fective trap frequency, reduced by the centrifugal "potential". 

In the thermodynamic limit (3Ml± tZ <C 1 and to = 0, 
G(x,0,0) = and (for a narrow FR, 7 <C 1) we re- 
cover the BCS-BEC crossover with T c (5) ranging from an 
exponentially small BCS-regime (5 > 2e F ^> fc^T) value 



- T BCS 



-2+c 



e F exp 



S-2ei 



■yep 



through T c c 



T C BEC [1 - ] 1/3 in the crossover < S < 2e F regime 

and saturating at T C BEC = (2C(3)y 1/3 hUN 1 ^ 3 for large neg- 
ative S. 

At finite u> and low T, /30j_ ^S> 1 (specializing for sim- 
plicity to an anisotropic trap f3fl± 3> 1, /3Q Z <§; 1), using 

G(x > l,a»l,k<l)wl - ic+iln(^),wefindan 



implicit equation /tyo 



c2 



2 e 



U-V2 

F 



quoted in the introduction. In the BCS and crossover regimes 
the implicit ui C 2 dependence enters through /.i as the solu- 
tion of the number equation, Eq. ( fl4l , that gives fi w € F « 
(SiV) 1 / 3 ^^]. The resulting solution for w C 2 is illustrated 
in Fig. [U with asymptotics, controlled by a dimensionless pa- 
rameter 3 ln(e F /M7jJ, summarized in Eqs.{T][2). 

We find that ui C 2 is driven to zero for S > S c m 2e F + 
^e F \n(e F /hil±). This corresponds to a vanishing of T" =0 
when the BCS condensation energy A 2 /e F becomes compa- 
rable to trap level spacing hil±. Equivalently, this condition 
corresponds to an oscillator length a± = ^/h/mfl± dropping 
to the coherence length £ = hv F /irA. 

In the opposite limit of low detuning, oj C 2(S -C S c ) — ► flj, 
where the centrifugal and trapping potentials nearly cancel, 
the system becomes translationally invariant and we obtain 
the bulk result hui C 2 ~ 8~ 1 e 2+c A 2 c2 /e'^ c2 , corresponding 
to Gorkov's (fixed chemical potential) prediction for type- 
II superconductors [8[|. This matches onto the BEC regime 
u C 2(5 < 0) ~ where a purely molecular superfluid can 
only be destroyed by a sufficient dilution down to comparable 
boson and vortex densities, driving a quantum transition into 
a variety of bosonic quantum Hall states |jl lj| . 

At high T (near T°) and slow rotation (ui <C w^ =0 ) in 
the BCS regime we expand G(x ^> l,o, 0) in the limit 
ax = (3 2 /ihn±(l + Ld 2 /Vt 2 ^) <C 1 (remaining in the de- 
generate limit, /3/i = x 1), obtaining a quadratic sup- 
pression of T c with rotation and trap frequencies given in 
Eq. (3). In the homogeneous limit to — > f2J, this reduces 

to T c M/T B b c u s lk « 1 - WftMAi/(* fl Tf <*) a , giving a lin- 
ear reduction with |u>|, as expected from Abrikosov's theory 
of H c2 (T) ft 

In the crossover regime at slow rotation, /i m S/2 < e F 
and a finite fraction of the atomic Fermi sea is bound into 
a molecular superfluid, giving N f (uj,S) w (5/2e F ) 3 N = 
(l/3)(<5/2MT) 3 /[l - (w/Ol) 2 ] < AT. Because_ the atomic 
density of states increases with according to 1 /S7[w] 3 (head- 
ing toward the extensive degeneracy of Landau levels in the 
u) — > limit), at sufficiently large w, e'p. drops below 5/2. 
The crossover frequency w*(<5), defined by e^* = (5/2, corre- 



sponds (for a narrow resonance 7 < 1) to a vanishing molec- 
ular condensate, N F (i^*,S) — N and is given by Eq. ©. 
w*(<$) is a lower bound for lu C 2(5) and marks the crossover 
between a strongly-paired molecular superfluid, with T c w 

T c BEC [JV B (a;,J)] 



r e (w,*)«T, 



BEC 



(5 
27^ 



n 1/3 



or 



, for w < u>*(6), 
(15) 

and a weakly-paired BCS superfluid with T c w T BCS (e£) 
for w*(5) < a; < U! C 2(S). For 5 < all atoms bind into 
molecules, N F (/1 < 0) is exponentially small, and the above 
result crosses over to the BEC value, T c (u>, 0), Eq. ©, quoted 
in the introduction. The crossover at uo w w»(<5) (or equiva- 
lently at S w 2e^*) can be seen in T c (<5, w) displayed for a full 
range of detunings and rotation rates in Fig. [2] 

To summarize, we have studied a resonantly-interacting 
atomic Fermi gas across a FR and computed uj C 2 (8, T) below 
which the gas becomes unstable to a rotating paired super- 
fluid. 
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